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We examine numerically vortex families near band edges of the Bloch wave spectrum in 
the Gross-Pitaevskii equation with a two-dimensional periodic potential and in the discrete 
nonlinear Schrodinger equation. We show that besides vortex families that terminate at a 
small distance from the band edges via fold bifurcations there exist vortex families that are 
continued all way to the band edges. 

I. INTRODUCTION 
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Many physical problems in periodic media with the Kerr (cubic) nonlinearity are governed by 
the Gross-Pitaevskii equation with a periodic potential. Examples are Bose-Einstein condensates 



n 



ized states in 



23| 



in optical lattices 13] and photonic-crystal fibers 18(. Interest in properties of local: 
this model has stimulated a number of mathematical works devoted to this subject 

An interesting problem that arises in this context is the possibility of bifurcations of stationary 
localized states from edges of Bloch bands in the wave spectrum of the Schrodinger operator with 



a periodic potential 



First pioneering works in this direction were completed by C. Stuart and 



his students [3, Q, 2l|. In the physics literature the asymptotic approximations of gap solitons 



bifurcating from band edges were developed by various authors in one dimension 
dimensions ^, Q, [2^ . 



id ] and two 



m 



In one dimension it was discovered numerically in [23|, Chapter 6.2] and explained analytically 
Q] that while single-pulse gap solitons bifurcate continuously from band edges, double-pulse gap 



solitons do not bifurcate from the band edges but experience fold bifurcations at a small distance 
from the edges. The situation becomes even more interesting in the space of two dimensions, where 
besides gap solitons, vortex solutions are possible in periodic potentials Q, I^, 



contradiction arises between the analytical results of 



221 ]. However, a 



SJ suggesting a continuous family of the 



fundamental vortex solutions bifurcating from the band edges and the numerical results of [23] (see 



2 

also [23', Chapter 6.5]) suggesting a fold bifurcation of vortex families at a small distance from the 
band edges. This contradiction will be inspected in this paper by using numerical computations. 

We will show that there do exist continuous families of the fundamental vortex solutions bi- 
furcating from band edges, according to the theory in 0, Numerical approximations of these 
families near band edges suffer, however, from the fact that the vortex localization is too broad and 
hence extends beyond the chosen computational domain. As a result, a spurious fold bifurcation 
occurs for the fundamental vortex family before the family reaches the band edge. If the size of 
the computational domain is enlarged, the location of the spurious fold bifurcation moves closer to 
the band edge. At the same time there are other vortex families, found also in 2^, which feature 



a true fold bifurcation at a small distance from a band edge. The fold location is independent of 
the size of the computational domain for these vortex families. 

The paper is organized as follows. Section |TT] introduces the two models which we in- 
spect, namely the Gross-Pitaevskii equation with a periodic potential and the discrete nonlinear 
Schrodinger equation. The connection between these models as well as the asymptotics of gap 
solitons near the band edges are reviewed. Section IIIII gives numerical results for the family of 
fundamental vortices. Section [TVl illustrates fold bifurcations for families of quadrupole and dipole 
vortex configurations. In Section |V] we summarize our findings. 



II. MODELS 



The stationary Gross-Pitaevskii equation with a periodic potential in the space of two dimen- 
sions takes the form 

-A^ + "\/(x,2/)(^-|(^|V = ^'^, (x,y)EM^ (II.l) 

where the focusing case is considered, the 27r-periodic potential V in each coordinate is assumed 
to be bounded, and G M is taken in a spectral gap of the Schrodinger operator 

Lo := -A + V. 

For simplicity, we assume that the periodic potential V has even symmetries with respect to reflec- 
tions about X = and y = 0. 

When u is close to the upper edge ujq of a spectral gap of Lq, a slowly varying envelope 
approximation of localized states can be derived and rigorously justified in the focusing stationary 



Gross-Pitaevskii equation . In the simplest case when uq is attained by only one extremum 
of the band structure and the Hessian at the extremum is definite, the resulting approximation is 

(p{x,y) = e'il;{ex,ey)(po{x,y) + OH''{e^^'^), w = wq - e^, (II. 2) 

where s > 1 is arbitrary, ifQ is the Bloch function at the band edge ojq, and ^ = il'{X,Y) in slow 
variables X = ex and Y = ey satisfies the stationary nonlinear Schrodinger (NLS) equation. This 
effective NLS equation is written in the form, 

ai^xx + i'YY)+P\i'\^i^ = i^, (II.3) 

where a > is related to the band curvature at the point ojq and /3 > is related to a norm 
of the Bloch function ipQ. We note that the leading-order term eip{ex,ey)ipo{x,y) has the order 
(Dh''{^) as e — )■ and hence expansion pi.2p shows that the perturbation term is smaller than the 
leading-order term in the norm, where is Sobolev space of square integrable functions and 
their derivatives up to the s-th order. When loq is the (upper) edge of the semiinfinite gap, the 
error was shown to be Oj/ife) or OL°°{e'^) [9]. 

The main theorem of [5|, l6|] states that if ip satisfies certain reversibility symmetries such as 

V(X, Y) = ±i;{-X, Y) = ±V^(X, -Y) (11.4) 

or 

ij{X,Y) = ±i;{Y, X) = ±i,{-Y, -X) , (II.5) 

and if the linearization of the stationary NLS equation pi.Sp is non-degenerate, then a localized 
solution of the Gross-Pitaevskii equation pi.ip with the asymptotic expansion pi.2p exists in 
for this ip. In particular, the stationary NLS equation pi.Sp admits the fundamental vortex of 
charge m S N, 



V'(X,y) = p(i^)e''"^ R = yJX^ + e = arg(X + iy), (IL6) 

where p{R) > for all ii > satisfies a certain differential equation that follows from the stationary 
NLS equation pi.3p . Vortex solution (jTLBj) satisfies symmetry pi.4p and the linearization of (jTLSj) 
at this Tp is non-degenerate. Hence conditions of the main theorem in P, are validated and there 
exists a unique localized solutions of pi.ip continued from this ip with the asymptotic expansion 
pi.2p . Continuation of fundamental vortices is considered in Section IIIII 



In the tight-binding Hmit of narrow spectral bands, several authors . [isl rigorously justified 



that the localized states of the stationary Gross-Pitaevskii equation (|II.ip can be described by the 
localized states of the stationary discrete nonlinear Schrodinger (DNLS) equation, 

- (Adisc0)m,n - \4'm,n\'^4'm,n = (^4'm,n, {m,n) G Z^, (II. 7) 

where 

(^disc0) m,n — 4'm+l,n m,n+l ~l~ (Pm—l,n ~l~ (Pm,n—l ^4'm,n 

and (jj ^ (t(— Ajisc) = [0,4]. The DNLS equation can simplify numerical approximations of the 
continuous Gross-Pitaevskii equation but does not change properties of localized states. In par- 
ticular, bifurcations of localized states are possible from the band edge w = in the focusing case. 
Moreover, the same method of asymptotic multi-scale expansions can be adopted to the DNLS 
equation with the expansion 

0m,n = eilj{em,en) + op(l), uj = -e^, (II-8) 

where "0 satisfies the same stationary NLS equation pi.3p . A rigorous justification of the contin- 
uous NLS equation as an asymptotic model for ground states of the DNLS equation was recently 
developed in Q]. Approximations for vortices in this context were not obtained to the best of 
our knowledge. 

III. VORTEX FAMILY CONNECTED TO THE SPECTRAL EDGE 

We compute here a family of the fundamental vortices in the DNLS equation pi.7p by using 
the near-edge asymptotics pi.SP with ^ given by the continuous vortex (|II.6P . We will also compare 
this behavior with the one in the Gross-Pitaevskii equation pi.ip . 

Choosing the vortex of charge one in the form pi.6P with m = 1, we compute the positive 
spatial profile p by the shooting method. We let next e = \/0.03, so that the expansion (jll.SP 
produces an initial guess for a solution (p of the DNLS equation 01. 7p with uj = —0.03 and compute 
(j) via the Newton's method. 

Next, we continue the family in the (co, ||(/>||p)-plane using the pseudo-arclength continuation 



Id, 



111 ], in which both and uj are unknowns, combined with the Newton's method. The resulting 
solution family is plotted in Figure (H^a) using the computational domain [—42,42]^ C Z^. The 
starting point at uj = —e^ = —0.03 is marked as D in Figure [^b). 
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The family of vortices with charge m = 1 seems to fold and never reach the spectral edge 
contrary to the approximation pi.Sp . In Figure [T]^b) this folding is, however, shown to be merely a 
numerical artifact caused by the truncation of the infinite domain Z^. The fold location approaches 
the edge uj = ojq = as the computational domain is enlarged. The family branch containing A — D 
thus terminates at the edge if computed on domains of diverging size. 

(a) 
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FIG. 1: Family of vortex solutions of (|II.7|) continued from the vortex (jll.6[) via the envelope approximation 
(jIL8| at a; = —0.03 (point D). (a) A fixed computational domain is used with = 85. (b) Detail of the 
vicinity of the spectral edge for a range of sizes of the computational domain. 

In the vicinity of the fold bifurcation, the solutions are inherent to the truncated domain and 
do not correspond to any solution of the DNLS equation pi.7p on Z^. The fold divides the family 
into two branches. Only the first branch containing points A,B,C,D is a continuation of the 
vortex family with the asymptotics given by pi.SP with ip as the continuous vortex (jTLGj). The 
other branch corresponds to a different family of vortex solutions. 

Figure [2] shows the solutions labeled ^ — G in Figure [H The limiting behavior for u — )• — oo 
along the branch with A and i? is a vortex with a square structure with sides of length five and 
three active sites on each side, i.e. the excited sites are at 

(-1, -2), (0, -2), (1, -2), (-1, 2), (0, 2), (1, 2), (-2, -1), (-2, 0), (-2, 1), (2, -1), (2, 0), (2, 1) e 

Along the other branch beyond the point G in the direction of decreasing uj the four solution peaks 
get further localized approaching single site excitations as w — )• —oo. The vortices keep charge one 
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FIG. 2: A-G: Modulus of the discrete vortex solutions labeled in Fig. [T] Bottom right: plots of the complex 
phase for vortices A and G. 

along the whole family. The complex phase for vortices A and G is plotted on the bottom right 
panel of Figure [2 

A similar situation arises for vortex solutions of the continuous Gross-Pitaevskii equation . 
In Figures [3] and [U we present a family of vortices for 

V{x,y) = 6sin^(x) + 6sin^(y), 



which was the potential used in [22]. We consider the vicinity of the lowest spectral edge coq w 
4.1264. The selected family is qualitatively similar to the discrete vortex family above and also 
terminates at the gap edge. The computational domain is [—L/2, L/2]'^ C and the stationary 
equation pi.l|l is discretized via central difference formulas of order 4. 



22( 1 that no vortex families can be continued 



The results of this section contradict the claim from 
to the 
from 

symmetries 01. 4^ or (jll.Sp and the non-degeneracy conditions, e.g., for the fundamental vortex 
solutions. 



band edge of the Bloch spectrum. These results illustrate the validity of the main theorems 
[s,^], which point out the possibility of such continuations for solutions satisfying reversibility 
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FIG. 3: Family of vortex solutions of (III.1[) continued from the vortex (jll.6[) via the envelope approximation 
pi.Sp at cj « 4.09 (point B). (a) A fixed computational domain is used with L = 93. (b) Detail of the vicinity 
of the spectral edge for a range of sizes of the computational domain. 

IV. QUADRUPOLE AND DIPOLE VORTEX FAMILIES DISCONNECTED FROM THE 

SPECTRAL EDGE 

A classical example of a vortex solution of the Gross-Pitaevskii equation pi.ip and of the 
DNLS equation pi.7p is the quadrupole vortex with the four nearest neighbor sites excited in a 
square arrangement in the asymptotics oj — t- — oo, i.e. with the excited sites at 

(-1,0), (1,0), (0,-1), (0,1) GZ2. 

For the Gross-Pitaevskii equation the excited sites are understood as locations of the single wells 
(minima) of the periodic potential V. We obtain this family for the DNLS equation via the initial 
guess 

at w = —1.2 (point D in Figure [5]), where 5m,miSn,ni is the Kronecker unit vector at (mi,ni) on 
Z'^. The family is then continued, once again, via the pseudo-arclength continuation combined with 
the Newton iteration. 

In Figures [5] and El the family curve and several solution profiles are plotted. The computation 
was performed on the domains [— A^, A^]^ with N = 10,16,20,40 and the curves in the (cj, H^Hp) 
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FIG. 4: Modulus of the continuous vortex solutions labeled in Fig. [S] The inset in A shows a qualitative 
plot of the complex phase for all vortices A-F. 




FIG. 5: Family of discrete vortex solutions continued from a quadrupole vortex and four nearest neighbor 
excited sites as w — >■ — oo. 

plane remained within the distance of 0(10^®) for all these and did not approach the edge. This 
family folds and does not reach the spectral edge and all its solutions remain tightly localized. It 
does not, therefore, contain a slowly varying solution near the band edge, which could be approx- 
imated via the envelope approximation pi.8p . Qualitatively this family corresponds to the family 
o — d in Fig. 1(a) of 2^. The branch past the fold, i.e. the branch with A and B, has as the 



asymptotic profile for lo — )• — oo a square vortex with excited sites at 

(-1, -1), (-1, 0), (-1, 1), (0, 1), (I, I), (1,0), (1, -I), (0, -1) G 
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FIG. 6: Modulus of the discrete quadrupole vortex solutions labeled in Fig. [5j The inset in A shows a 
qualitative plot of the complex phase. 

We also consider a family of real dipole solutions (p which are odd in the n-index and satisfy 
0m, = for all m G Z. It is constructed via the "hand made" initial guess 

0m,n = l-5(5m,o5n,l - 1 .5(5m,o5n -1 , (m, n) € Z^, 

at a; = —1.2 (point D in Figure[7]). The family is plotted in Figures [7] and [HI The computation was 

2 N 



performed on the domains [— A^, A^]^ € with = 10, 15,20,40 and the curves in the (w, ||</'||?2) 



plane remained within the distance of 0(10^^) from each other for all these and did not approach 
the edge. This family contains again only tightly localized solitons and does not continue to the 
spectral edge either. Along the branch with points E and D the asymptotic profile for uj — t- — oo is 
the dipole with the excited sites (0, —1) and (0, 1). Along the other branch the asymptotic profile 
is twice as broad with the excited sites at 

(0,-2), (0,-1), (0,1), (0,2) ez2. 



22l |. where all quadrupole 



The results of this section confirm the previous numerical results from 
vortex families have a fold bifurcation at a small distance from the band edge of the Bloch spectrum. 
Solutions throughout our quadrupole and dipole families remain tightly localized so that they are 
not compatible with the slowly varying approximation pi.sp . It is, however, possible that there 
exist other solution families, which have a dipole or quadrupole asymptotic form as a; — t- — oo and 
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FIG. 7: Family of discrete vortex solutions continued from a dipole vortex with two nearest neighbor excited 
sites as w — >■ — oo. 
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FIG. 8: Modulus of the discrete dipole vortex solutions labeled in Fig. [71 



which do bifurcate from the spectral edge. The bifurcation would be guaranteed by the results of 
[5,^] if the continuous NLS equation ()II.3P had dipole or quadrupole solutions, which we are not 
aware of. 



V. CONCLUSION 



We have numerically demonstrated in both the Gross-Pitaevskii equation with a two- 
dimensional periodic potential and the discrete nonlinear Schrodinger (DNLS) equation that there 
are families of fundamental vortices bifurcating from spectral edges of the Bloch wave spectrum. 



This is in agreement with the analysis in and in contradiction with the claim in [22] that 

no vortex families continue to a spectral edge. Our fundamental vortex families complement the 



22 



23[ which all terminate at a distance from the spectral edge via 



vortex families constructed in 
a fold bifurcation. 

We have also investigated families of quadrupole and dipole vortex configurations of the DNLS 
equation. The selected families do terminate via fold bifurcations, which are located at a small 
distance from the band edge, independently of the size of the computational domain. It is an open 
question whether there exist other families of quadrupole and dipole vortex configurations, which 
bifurcate from spectral edges of the Bloch wave spectrum. 
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